Detecting ill posed boundary conditions in General Relativity 
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A persistent challenge in numerical relativity is the correct specification of boundary conditions. 
In this work we consider a many parameter family of symmetric hyperbolic initial-boundary value 
formulations for the linearized Einstein equations and analyze its well posedness using the Laplace- 
Fourier technique. By using this technique ill posed modes can be detected and thus a necessary 
condition for well posedness is provided. We focus on the following types of boundary conditions: i) 
Boundary conditions that have been shown to preserve the constraints, ii) boundary conditions that 
result from setting the ingoing constraint characteristic fields to zero and iii) boundary conditions 
that result from considering the projection of Einstein's equations along the normal to the boundary 
surface. While we show that in case i) there are no ill posed modes, our analysis reveals that, unless 
the parameters in the formulation are chosen with care, there exist ill posed constraint violating 
modes in the remaining cases. 
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' Obtaining a long time convergent numerical simulation of a binary black hole spacetime in domains with artificial 
boundaries continues to be a challenge in numerical relativity and one which has recently received a substantial amount 
of attention, notably in the case of hyperbolic formulations (see [1, 2] for reviews). The challenge remains in part 
| because of the difficulty in specifying boundary conditions. It has been recognized [3, 4, 5, 6, 7, 8, 9, 10, 11, 12] that 
the boundary conditions have to satisfy two important requirements. First, they have to preserve the constraints. 
By this we mean that they must guarantee that if the constraints are satisfied initially they are also satisfied at 
later times. We refer to boundary conditions that satisfy this property as constraint preserving boundary conditions 
(CPBC). Second, the boundary conditions have to be such that the resulting initial-boundary value problem (IBVP) 
is well posed. This means that given initial and boundary data a unique solution exists and that at each fixed time the 
solution depends continuously on the data. Well posedness is a necessary condition for the construction of consistent 
f~>) I and stable finite difference schemes [13, 14]. 

When the evolution equations are in symmetric hyperbolic form one usually specifics maximal dissipative boundary 
conditions [15]. Under certain technical assumptions, these conditions guarantee that the resulting IBVP is well posed 
[16, 17]. Using maximal dissipative boundary conditions, Friedrich and Nagy [3] were able to find well posed CPBC 
for a particular formulation of the full nonlinear vacuum equations. However, most of the hyperbolic formulations 
used in numerical relativity are based on evolution equations that use a different set of variables than in Ref. [3]. 
For these formulations, the derivation of well posed CPBC seems to be more difficult. Part of the problem stems 
from the fact that CPBC result in a set of partial differential equations that must hold at the boundary surface, and 
it is not always possible to cast these equations into the form of maximal dissipative boundary conditions. This is 
probably the reason why current well posed CPBC for formulations other than that used in Ref. [3] are either limited 
to homogeneous boundary data [9] or to linearizations around a Minkowski background [7, 8, 9]. Even in those cases, 
the CPBC obtained so far might be too restrictive in the sense that they do not allow the specification of the physical 
quantities at the boundary with the freedom one would like to have. For example, the well posed boundary conditions 
obtained in Ref. [7] involve a coupling between the in- and outgoing variables and, likely, this coupling will introduce 
reflections at the boundary. Therefore, more general techniques are desirable in order to show well posedness for more 
generic CPBC. 

In this article, we use the Laplace-Fourier technique to analyze boundary conditions in linearized General Relativity. 
This technique is very useful when the evolution equations are linear and have constant coefficients since it can be 
applied to boundary conditions that are more general than the maximal dissipative ones. Specifically, it can be applied 
to boundary conditions which have the form of differential equations at the boundary. Furthermore, the method is 
capable of detecting the presence of ill posed modes analytically. Ill posed modes are solutions to the IBVP that grow 
exponentially in time with an exponential factor that can be arbitrarily large, and their existence makes it impossible 
for the solution to depend continuously on the data. The Laplace-Fourier technique therefore provides us with a 
necessary condition for well posedness. However, it should be emphasized that the absence of ill posed modes (as 
defined in this article) does not automatically guarantee well posedness. Although more complicated in this case, 
results for the variable coefficient case are available by freezing the coefficients at the boundary (see [18, 19]). 

This article is organized as follows: The conditions under which the specification of non-maximal dissipative bound- 
ary conditions for symmetric hyperbolic systems with constant coefficients yields ill posed modes are reviewed in section 
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II. In section III we discuss the boundary conditions that have been considered for the generalized Einstcin-Christoffel 
formulation of Einstein's equations [20, 21] when linearized around flat spacetime. The generalized Einstcin-Christoffel 
system is a family of symmetric hyperbolic formulations that is parametrized by a constant rj. The boundary condi- 
tions we are considering are: i) The CPBC that were considered in Ref. [7] and that are based on solving a closed 
evolution system at the boundary and on maximal dissipative boundary conditions, ii) Boundary conditions that are 
obtained by setting the ingoing constraint characteristic fields to zero, iii) Boundary conditions that are obtained by 
considering the projection of Einstein's equations along the normal to the boundary surface, as recently proposed by 
Frittelli and Gomez [10, 11]. In section IV we apply the techniques discussed in section II and show that the cases 
ii) and iii) suffer from the presence of ill posed modes unless the parameter rj in the generalized Einstcin-Christoffel 
formulation lies in a specific range. We also show that there are no ill posed modes in case i) which is consistent with 
the well posedness estimates derived in Ref. [7] . In section V we show that the ill posed modes we have found in cases 
ii) and iii) do, in fact, all violate the constraints. This means that the evolution system for the constraint variables is 
ill posed in those cases. Since this system is always strongly hyperbolic and since our boundary conditions are con- 
structed from specifying maximal dissipative conditions for this system, this also illustrates that maximal dissipative 
boundary conditions do not necessarily yield a well posed formulation if the evolution equations are strongly hyper- 
bolic (but not symmetrizable) . In particular, our calculations show that the boundary conditions that are constructed 
following the schemes ii) and iii) do not necessarily lead to CPBC and that one should always check the evolution 
system for the constraint variables. Our results and their implications on deriving well posed CPBC are discussed in 
section VI. A similar analysis for the Frittelli- Reula [22] system has been undertaken by Stewart [4]. 



II. DETECTING ILL POSED MODES 



In this section, we review the techniques that can reveal the presence of ill posed modes. They are based on a 
Laplace transformation in time and on a Fourier transformation in the spatial directions that are tangential to the 
boundaries and are described in more detail in Refs. [13, 23]. For simplicity, we restrict the following discussion to 
the 2D case; the generalization to 3D is straightforward. 

Consider a 2D first order in time and space linear evolution equation of the form 

d t u = Ad x u + Bd y u 7 (1) 

where u = u(t,x,y) is a vector-valued function and the matrices A and B are constant and symmetric. We consider 
solutions to Eq. (1) on the domain t > 0, x > 0, — ir < y < +ir with initial data 

u(0,x,y) = f{x,y) (2) 

and boundary conditions at x = of the form 

L(d u d y )u{t,0,y)=g(t,y), (3) 

where L is a linear operator with constant coefficients that only involves derivatives which are tangential to the 
boundary. For technical reasons, we assume that L(d t , d y ) is homogeneous in the sense that L((idt, ftd y ) = /jL(d t ,d y ) 
for all positive fi. We also assume periodic boundary conditions in the y-direction (similar conclusions hold for the 
case — oo < y < +oo). 

The IBVP (1), (2), (3) is said to be well posed [24], if given smooth square intcgrablc data /, g there exists a unique 
smooth solution. Furthermore, there are constants C, a such that 



\\u(t,.)\\ 2 <Ce" 
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+ \\9(T,.)VdT 



(4) 



for alii > and all square intcgrable data / and g. Here, \\u(t, .)|| denotes the L 2 norm of u defined by ||it(t, .)|| 2 = 
J x>Q \u(t,x,y)\ 2 dxdy and similarly, ||/|| 2 = J x>Q \f(x, y)\ 2 dxdy and ||g(r,.)|| 2 = / \g(r, y)\ 2 dy. The estimate (4) 
implies that for each fixed t, the solution depends continuously on the data / and g. 

A first step in checking if a given initial-boundary formulation satisfies a well posedness inequality of the type (4) 
is to look for solutions of the problem with homogeneous data (g = 0) which are of the form 



u(t, Xl y) = e st+ ^yu(x), (5) 

where lu is an integer, s is complex with Re(s) > and u(x) is a smooth function that lies in L 2 (0,oo). If such a 
solution exists, the problem cannot be well posed. In order to see this we notice that the functions 

u m (t,x,y) = e m ^ st+i ^u(mx), (6) 
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where m = 1,2,3,... can be arbitrarily large are also solutions and since ||w m (£, .)||/||ti TO (0, .)|| — exp(mRe(s)i) the 
estimate (4) cannot hold with constants C and a that are independent of the initial data. Therefore, an obvious check 
for well posedness is to see whether or not Eqs. (1), (3) admit nontrivial solutions of the form (5) with homogeneous 
boundary data. 

Using expression (5) in Eqs. (1), (3), we obtain (for 5 = 0) 

su = Ad x u + itoBu, (7) 
L(s,icu)u(0) = 0. (8) 

These equations form a system of ordinary differential equations and can be solved analytically. In order to do so, we 
first bring A to diagonal form through an orthonormal transformation. Thus, the matrix B is still symmetric and we 
can write 



A = 




A x 



where A\ = diag(Ai, A p , A p+ i, \ p + q ) with Ai, ...A p real and negative and A p +i, X p + q real and positive. Here, 
p and q are the number of in- and outgoing modes, respectively. Accordingly, we write 

D / -Boo B i \ , . / u 

B= {b 10 B 11 )> L = ( L o,Li), u ={u 1 

Now the zero components of Eq. (7) yield the following algebraic relation between u and u\: 

Soouo = —Soiui , (9) 

where we have introduced the matrix S = S(s,cu) = si — iu)B. Since the matrix B is symmetric, the matrix 
Soo = si — iu)B o is invertible for all Re(s) > and all integer ui, and we can express uo in terms of 

u = -SqqSoiUi ■ (10) 

Inserting this into the remaining equations of the system (7), (8), we obtain the reduced problem 

d x ui = M(s,w)u! , (11) 
Lui = 0, (12) 

where 

M{s,uj) = A7 1 (Su(s,u) - SwS^Soiis,^)) , 
L(s,w) = Li—LqSqqSqi- 

One can show [19] that for Re(s) > the matrix M(s,uj) has exactly p eigenvalues with negative real parts and 
exactly q eigenvalues with positive real parts (the eigenvalues are counted according to their algebraic multiplicity.) 

The eigenvalues of M that have positive real part lead to exponential growth in x. Since the solution u has to be 
in L 2 (0, oo), the integration constants have to be chosen such that there is no such growth in x. In order to achieve 
this, we choose, for each (s,uj), a unitary matrix U = U(s,ui) that brings M{s,w) in upper triangular form: 



U(s,u) 1 M(s,w)U{s,uj) =y V Q ' > M ° V (s ' )W ) J • 



(13) 



Here, M_ (M + ) is an upper triangular matrix whose eigenvalues have negative (positive) real parts. If we introduce 
the new variable v(x) = U(s,lj)~ 1 ui(x), system (11), (12) becomes 

d x v^(x) = M-(s,lj)v-(x) + M {s,u>)v + (x) , 
d x v + (x) = M + (s,uj)v + (x) , 
i_w_(0)+i+w+(0) = 0, 

where (£-, L + ) = LU. It follows that v+ must vanish for v to be in £ 2 (0, oo). This implies that V-(x) = exp(M_x)a_ 
where cr_ has to satisfy the boundary condition i_cr_ = 0. We conclude that the system (7,8) has only the trivial 
solution if and only if the determinant condition[25] 

det£_(s,w) ^ 0, Re(s) > 0, (14) 
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is satisfied. (In particular, _L_ must be a square matrix of dimension p. This means that we need exactly as many 
independent boundary conditions as there are ingoing modes) . If the determinant condition is violated at some point 
(s,u>) — (sq,u>q), it is also violated for (s,u>) = m(so,ujo) with m — 1,2,3,... and the initial-boundary formulation 
admits solutions of the form (5) that grow exponentially in time where the exponential factor s can have arbitrarily 
large real part. 

In section IV we will discuss the determinant condition for the case of the linearized Einstein equations with 
boundaries. 



III. BOUNDARY CONDITIONS FOR THE LINEARIZED EINSTEIN-CHRISTOFFEL SYSTEM 

In this section we discuss boundary conditions for a linearization of the generalized Einstein-Christoffcl vacuum 
equations [20]. This formulation has the attractive feature that when linearized around flat spacetime written in 
Minkowski coordinates it simply reduces to a set of six wave equations, written in first order form: 

dtKij = -S kl d k fu 3 , (15) 
d t fk t3 = -d k Kij . (16) 

Here, K i3 denotes the linearized extrinsic curvature and the symbols fkij represent linear combinations of the linearized 
Christoffcl symbols Fkij'. 

fkij = F(ij)h + (^ki^[ SJ ]r + SkjT[si]r + ^ ^ %r[ sfc]r ^ . (17) 

The value of rj (which has to be different from zero) parametrizes the family of formulations. The particular case with 
j] = 4 corresponds to the original Einstein-Christoffel system derived by Anderson and York [21]. We set the shift 
to zero, and the lapse is linearized in such a way that it satisfies the densitized lapse gauge condition a = ^fg up to 
second order corrections, where g denotes the determinant of the three metric. A solution to the system (15), (16) is 
a solution to the linearized Einstein equations if and only if the constraints are satisfied. In terms of the constraint 
variables 

C = 1 S rs d r v s , (18) 
Cj = S rs {d r K sj -djK rs ), (19) 

77 — 4 

Ci kij = 2d {l f m +vd[iS k]{l v 3) + Sijdyvq , (20) 

where Vk — S^ifkij — fijk), the constraints are given by C — 0, Cj = 0, Cikij = 0. 

We consider the evolution system (15), (16) on the domain t > 0, x > 0, — n < y,z,< +n and introduce the 
characteristic variables in the x dircction[26] 

4 7) = +/*«)> (21) 

4 +) = ^(Kij-fxa), (22) 

«S = fMj ■ (23) 

Here and in the following capital Latin indices stand for the tangential directions y and z. When written in terms of 
these variables the evolution equations (15), (16) take the form 

= -d x utf ^S AB d A u% , (24) 



V2 
V2 

d t u% = -j=0 A ^+^), (26) 



= +d x u^ ^ AB dAU% , (25) 



and we see that the matrix A in Eq. (1) is diagonal. 
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For the constraints to be satisfied everywhere, when boundaries are present, one has to ensure that they are satisfied 
initially and that no constraint violating mode enters the domain. In order to ensure this, we follow the analysis in Rcf. 
[7] and first consider the evolution of the constraint variables with respect to the flux defined by the main evolution 
equations (15), (16). One can show [7] that the traceless part of Ctkij is constant in time, while the remaining 
constraints propagate according to 



d t C = j5 l ^d l C , (27) 



d t Cj = ^p-djC -8 rs d r T sj , (28) 
d t Tij = -diCj + (l - djd + | 6 l3 6 rs d r C s , (29) 



dtVn 



f - 3 ) d H C i] > ( 3 °) 



where = S rs (C r ij S + CV, rs ), and = S rs Cij rs . Introducing k = 1 — 3rj/A and the variables 

Qj = T i:j + 2r] - SijC = 5 rs (d r fi js - djfirs) + KdjVi - n5i 8 rs d r v s , (31) 



the characteristic fields can be written as [27] 



= ^=(c j + c xj ), m 
vj +) = ±=( c -c,). m 



V A f = C Aj + K (6 xj C xA -S Aj C xx ), (34) 



= - : ^C [ij] + (« + l)V tJ . (35) 

If the system (27), (28), (29), (30) is symmetric (or symmctrizablc) hyperbolic, one can guarantee that if the constraints 
are satisfied initially and if homogeneous maximal dissipative boundary conditions are given, the constraints will be 
satisfied everywhere. Therefore, we consider boundary conditions at x = which are of the form 

Vj") - aVj+) - 0, v{- } - bV { / ] = 0, (36) 
U L ) - c«W = 9xx, u a ~b - du { AB = 9AB, (37) 

where the magnitudes of a, 6, c and d are smaller or equal to 1 and iS A g = u^ AB — ^5 A bS cd Uq B denotes the traceless 

part of u'ab- I n order to express the conditions (36) in terms of the main variables Kij, fkij, we use the definition of 
the constraint variables, Eqs. (18), (19), (31), and the evolution equations (15), (16) in order to trade x-derivatives 
by time and tangential derivatives: 
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±d t u™ + d A u™ ± ^ d c u^ AB ± \d A (V2S CD u^ DB V2S-u% + «&-,,', ;/ ^ 



V2 



It follows from the energy estimates derived in Ref. [7] that when < r\ < 2 the conditions (36) guarantee that the 
constraints are satisfied everywhere if they are satisfied initially. In the following, we will also consider other values 
of f] and show that one might have ill posed modes if the parameter rj lies outside the interval (0,2). Notice that 
the conditions (36) do not involve derivatives normal to the boundary (d x ). They can be interpreted as evolution 

equations for the variables 5 AB (u AB + au^g) and u^ xA + i>u^ A at the boundary. The functions g xx and g A B are data 
that can be given freely for a combination of the in- and outgoing gauge and physical variables, respectively [7]. 
In the next section, we will analyze the following choices of parameters: 

1. a = — 1,6=1, c = d = 1: 

This corresponds to the Neumann boundary conditions that we have discussed in Ref. [7]. In this case the 
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boundary conditions can be recast in a closed evolution system at the boundary. Its solutions provide boundary 
data for the main evolution system in the form of maximal dissipative boundary conditions. When < r\ < 2 
one can derive well posedness estimates for the resulting IB VP and the boundary conditions can indeed be called 
CPBC. 

2. a = 1, b = -1, c = d= -1: 

This corresponds to the Dirichlet conditions specified in Ref. [7] . They can also be recast in a closed evolution 
system at the boundary and for < rj < 2 one has a well-posed IB VP with CPBC. 

3. a = 0, b = 0: 

This corresponds to setting the ingoing constraint variables to zero and might be the most obvious choice for 
obtaining CPBC. However, we will show in the next section that the resulting IBVP possesses ill posed modes 
unless the parameter r\ is chosen appropriately. 

4. a = 0, b= 1: 

These are the conditions that one obtains after linearizing the boundary conditions that were recently proposed 
in Ref. [11]. There, the Einstein-Christoffel formulation (77 = 4) is considered and the boundary conditions 
(36) are obtained by projecting Einstein's equations along the normal to the boundary rather than by analyzing 
the evolution of the constraints. In fact, one can show that setting G x a to zero and rewriting[28] the resulting 
equations in terms of the variables and fkij is equivalent to the second equation in (36) with 6=1, while 
setting G xx (Gt x ) to zero is equivalent to the first equation in (36) with a = 1 (a = — 1). In Ref. [11], the 
authors propose to set the linear combination G xx — Gt x to zero which would correspond to using a = in 
(36) . In the next section, we show that the resulting boundary conditions yield an ill-posed formulation if the 
parameter rj is not chosen appropriately. 



IV. LAPLACE-FOURIER ANALYSIS 



Following the analysis described in section II, we look for solutions to Eqs. (15), (16), (36), (37) with homogeneous 
boundary data (g xx = 0, c/ab = 0) and which are of the form 

u\f\t,x,y,z) = e st+i "v y+iw * z v§\x), (38) 
u%(t,x,y,z) = e^y+^u%(x), (39) 

where s is a complex number with positive real part and uj v and uj z are integers. For the solution to be square 

integrable, we require the functions u^p (x) and u%(x) to be in L 2 (0,oo). From Eqs. (38), (39) and (26) we obtain 
an algebraic condition 

su%=-^ A (uV + ^) (40) 

that can be used to eliminate the variable u^ij from the remaining equations. Inserting Eqs. (38), (39), (40) into 
Eqs. (24), (25) yields the ordinary differential equation 



M{s,w) % , (41) 



where 





M(s,w) = \ 2 2, 2« 2 (42 ) 

V s + 17 / 

and to — (u> y ,uj z ). The matrix M(s,of) has the eigenvalues ±-\/s 2 + w 2 . 

We first look at the case u = 0, which corresponds to solutions that have trivial y and z dependence. For those, 

the matrix M(s,u) is diagonal and since Re(s) > we see that we must have = for the solution to be in L 2 . 
The boundary conditions (36), (37) yield 

C } (0)-0, fi!ri(0)=0, ^ B «^(0) = 0, ^(0) = 0, (43) 
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therefore we have only the trivial solution. There are no ill posed modes with trivial dependence on the variables 
that are tangential to the boundary. We show now that the situation becomes rather more complicated when one 
considers modes that have a nontrivial tangential dependence. 

Assume that w^O. Following the analysis of section II we introduce a unitary matrix U = U(s,lo) that brings 
the matrix M{s,u) into upper triangular form. To lighten the notation we introduce the quantities ( = s/\w\, 
^ = yl + C 2 ! V'(C) = — C) 2 and N = 1 + |^i>(£)| 2 . One can then verify that the matrix (a star denotes complex 
conjugation) 

U(0="- 1/2 [l^ C £ J (44) 



is unitary and satisfies 



no-.«no = u(- ^°). WO- 1 **"* 2 *.?**') . («) 



where in A the branch is chosen such that for Re(() > 0, Re(A) > 0. In terms of the new variables [v\- \v^) T — 
U- 1 (UiJ ) Mj ) ) T Eq- (41) yields 

|w|~ 1 ax«{,~ ) = -A<4r)+M (C>i + \ (46) 
Iwl" 1 ^ = A^. (47) 

For the solution to be in L 2 , we must have vj^ = 0. This implies that v\~j ' = e~ x \-\ x aij, where Oij are constants 

which describe the value that v\- ' takes at the boundary. Using the matrix U(() we can express the u variables at 
the boundary as 

^(0) = -N^-fa^ (48) 
< +) (0) = N-V^MQvv , (49) 



u 



c 

(0) m-, _ i u A A7 ._i /2 |CI 

Aij 



(o) = C V-M))°«> ( 5 °) 



where a>yi = 

Using Eqs. (38), (48) and (49) in the boundary condition (37), we find that 

(1 + cyj(())a xx =0, (1 + #(C))<xab = 0, (51) 

where uab denotes the tracefree part of (jab- Since the function V(C) maps the half plane Re(£) > to the interior 
of the unit circle and since \c\ < 1, \d\ < 1, it follows that a xx = and &ab = 0. 

Next, we insert all of this into the boundary conditions (36). The result is more conveniently expressed if one 
introduces a normalized vector £4 that is orthogonal to uja and considers the components a XL0 = S AB a x A^B and 
&xi = S ab <t x a£,b- The projection of the second equation in (36) along £ implies that a x ^ must vanish, while the 
remaining equations in (36) imply that a = 5 AB (tab and cr XLU must satisfy the following 2x2 system: 



T , n ( a \_ . ( 2\(1 + aiP) - K (i + a)(A-C) 2i(l + aip) + in{l + o)(l + ip) 

^>\ct xu J- U ' + -i«(l + 6)(l + V) 2A(1 + 6^) + 2k(1 + 6)(A-C) 



(52) 



The determinant of L-(C) is given by 

det L_(C) = (6 + 4C 2 ) [(1 + o^(C))(l + &1KO) - ^(1 + «)(! + &M0] ■ (53) 

Clearly, the first factor cannot be zero since Re(£) > 0. Therefore, detL_(C) can only vanish if the term inside the 
square brackets does. 

We now focus on the different cases discussed in the previous section: 
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1. a = — 1, 6=1, c = d = 1: 

In this case, the second term inside the square brackets vanishes and the first term is never zero since [VKOI < 1- 
Therefore, the resulting formulation possesses no ill posed modes. Of course, when < rj < 2 this is consistent 
with our calculation in Ref. [7] where the estimates we have derived exclude the presence of ill posed modes. 

2. a = 1, 6 = -1, c = d = -1: 

The result is the same as in the previous case. 

3. a = 0, b = 0: 

In this case, the terms inside the square brackets simplify to 1 — K 2 tp((). A small calculation reveals that this 
can only be zero if £ = (n 2 — 1)/2|k| and k ^ 0. Therefore, detL_(£) has a zero with Re(C) > if and only 
if k 2 > 1. This is equivalent to rj < or rj > 8/3. Therefore, setting the ingoing constraint variables to zero 
in the family of generalized Einstein-Christoffel systems does indeed yield ill posed boundary conditions if the 
parameter rj lies outside the interval [0, 8/3]. 

4. a = 0, b = 1: 

Here the terms inside the square brackets reduce to 1 + (1 — 2k 2 )i/>(£). Since the function ip maps the positive 
real axis onto the open interval (0, 1) it follows that this expression never vanishes if and only if n 2 < 1. In 
particular, one has ill posed modes when 77 = 4 and the boundary conditions that were proposed in Ref. [11] 
yield, at least when linearized around flat spacetime, an ill posed initial-boundary formulation. On the other 
hand, if at the boundary one considers the equations G xy — G xz = and the equation G x t = instead of the 
combination G xx — G xt = one has a = — 1 and the resulting formulation does not in fact suffer from possessing 
ill posed modes. 



V. VIOLATIONS OF THE CONSTRAINTS 

In this section, we show that the ill posed modes we have found in the previous section violate the constraints. In 
order to see this, we use these ill posed modes to compute the constraint variables Cj. From = (u[j ^ +u^)/v / 2, 
Eqs. (19), (48), (49) and a xx = 0, &ab = 0, we have 

-' CI (l-^(C))(Aa + 7a^)exp[|w|(Ct-Ax + « C 2; A x A )] , (54) 




— , (1 - V>(0) {io + 2Xa XUJ ) cxp [y ((t - \x + iu A x A )} , (55) 
oJS s 

where (a, a xul ) is a nontrivial solution to Eq. (52). Since \ip(()\ < 1 for Re(£) > 0, and since 

det ( i 2A ) = 3 + 2C - 2 * °' ( 56 ) 

we see that the variables C x and uj a Ca cannot simultaneously vanish. Therefore, all the ill posed modes we have 
found are constraint violating modes. This means that under generic small perturbations of the initial data these 
modes will be excited and the constraint variables will grow exponentially with an exponential factor that can be 
arbitrarily large. In this sense, the boundary conditions that lead to ill posed modes do not preserve the constraints. 
We point out that the constraint variables constructed from any solution of the main evolution system (15), (16) with 
boundary conditions (36), (37) provide a solution of the evolution of the constraint variables, Eqs. (27), (28), (29), 
(30) with boundary conditions (36). Since we have shown that the constraint variables constructed from ill posed 
modes are ill posed modes themselves (see Eqs. (54), (55)), the IBVP for the constraint variables cannot be well 
posed. This emphasizes the importance of looking at the evolution system for the constraints and checking its well 
posedness when deriving CPBC for Einstein's equations. 

We conclude this section with two remarks. First, one can check that the evolution system for the constraint 
variables, Eqs. (27), (28), (29), (30), is strongly hyperbolic for any nonvanishing value of the parameter 77. On the 
other hand, our analysis above shows the existence of ill posed modes when r\ lies outside of the interval [0,8/3] 
and the coupling constants a and 6 are chosen as in the cases 3. and 4. of the previous section. This illustrates 
that applying maximal dissipative boundary conditions to evolution systems that are strongly hyperbolic (but not 
symmetrizable) does not necessarily yield a well posed problem. 

The second remark concerns the choice a = b = — 1 for the coupling constants in Eq. (36). In this case it follows 
that the determinant condition is always satisfied, regardless of the value for the parameter rj. In fact, one can show 
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that the resulting boundary conditions are constraint preserving: The evolution equations imply that the constraint 
variable Cj satisfies the wave equation: 

d\Cj = S rs d r d s Cj . (57) 

On the other hand, the choice a = b = — 1 corresponds to imposing the momentum constraint at the boundary. Since 
Cj satisfies the wave equation, this implies that Cj = everywhere, if Cj is satisfied initially. It then follows from 
Eqs. (27), (29), (30) that the remaining constraints are also satisfied if they are satisfied initially. This explains why 
one has CPBC for all rj ^ when a = b = — 1. However, the above argumentation is expected to break down when 
one considers the nonlinear regime since in this case lower order terms might prevent one from obtaining a closed 
system for Cj alone. In this case, one has to rely on the symmetrizer for the system (27), (28), (29), (30) which was 
constructed in Ref. [7], and one might not be able to show that the constraints propagate when r\ lies outside the 
interval [0, 2), even when a = b = —1. 

VI. CONCLUSIONS 

We have analyzed ill posed modes in the family of the generalized Einstein-Christoffcl formulation of Einstein's 
equations with boundaries. We considered boundary conditions which result from coupling the ingoing characteristic 
constraint variables to the outgoing ones. Specifically, the cases we have studied include the boundary conditions we 
have obtained in Ref. [7] and the boundary conditions that originate from considering the projection of Einstein's 
equations along the normal to the boundary. When linear fluctuations around Minkowski space are considered, we 
have shown that the formulation is subject to constraint violating ill posed modes unless the parameters in the 
equations and the coupling between the in- and outgoing constraint variables are chosen carefully In fact, it is not 
difficult to show that if the coupling constants a and b are real and satisfy — 1 < a < 1 and — 1 < b < 1 there are always 
ill posed modes as long as the parameter r\ lies outside the interval [0,8/3]. In particular, this is the case when the 
ingoing constraint variables are set to zero. Furthermore, there are ill posed modes for the boundary conditions that 
were obtained in Ref. [11] when applied to the linearized Einstein-Christoffel system (rj = 4). However, our analysis 
also reveals that these ill posed modes could easily be avoided by imposing a different linear combination of Einstein's 
equations at the boundary or by changing the parameter r\ such that it lies in the interval < r\ < 8/3. In any case, 
our analysis highlights the importance of studying the evolution system for the constraint variables and ensuring its 
well posedness since all the ill posed modes we have found are constraint violating. In particular, the formulations we 
have studied in this article show that even though the main evolution system is symmetric hyperbolic, the evolution 
equations for the constraint variables is not necessarily symmetrizable. For the cases in which the propagation of 
the constraints is described by a system that is strongly hyperbolic (but not symmetrizable) we have shown that 
specifying maximal dissipative boundary conditions can lead to an ill posed system. 

It is interesting to note that all the ill posed modes that appear have a nontrivial dependence in the spatial directions 
that are tangential to the boundary surface. Therefore, such modes would not be present in the one-dimensional case. 
This might explain why the numerical simulations in Ref. [6], where the Einstein-Christoffel system (rj — 4) was 
evolved using boundary conditions obtained by setting the ingoing constraints to zero, did not show any ill posed 
modes. 

The simple analytic method we have used in this article, which is based on the determinant condition (14), should 
be used to test the well posedness of the boundary conditions before numerically evolving any evolution system since 
the presence of ill posed modes would detrimentally affect numerical stability. However, we also stress that more 
work is required to derive sufficient conditions for well posedness for the choices of parameters when the determinant 
condition is satisfied. In particular, it would be worthwhile to analyze CPBC where the incoming physical variables 
can be freely specified. 

Acknowledgments 

We particularly thank Gabriel Nagy for many valuable discussions. We also thank L. Lehner, J. Pullin, O. Reula 
and M. Tiglio for useful comments. This work was supported in part by grants NSF-PHY-9800973, NSF-INT-0204937, 
and the Horace C. Hearne Jr. Institute of Theoretical Physics. 



[1] O. Reula, Living Reviews in Relativity 1, 3 (1998). 



10 



[2] L. Lehner, Class. Quantum Grav. 18, R25 (2001). 

[3] H. Friedrich and G. Nagy, Comm. Math. Phys. 201, 619 (1999). 

[4] J.M. Stewart, Class. Quantum Grav. 15, 2865 (1998). 

[5] M.S. Iriondo and O.A. Reula, Phys. Rev. D 65, 044024 (2002). 

[6] G. Calabrese, L. Lehner, and M. Tiglio, Phys. Rev. D 65, 104031 (2002). 

[7] G. Calabrese, J. Pullin, O. Sarbach, M. Tiglio, and O. Reula, "Well posed constraint preserving boundary conditions for 

the linearized Einstein equations," arXiv:gr-qc/0209017. 
[8] B. Szilagyi, B. Schmidt, and J. Winicour, Phys. Rev. D 65, 064015 (2002). 

[9] B. Szilagyi and J. Winicour, "Well Posed Initial- Boundary Evolution in General Relativity," arXiv:gr-qc/0205044. 
[10] S. Frittelli and R. Gomez, "Einstein boundary conditions for the 3+1 Einstein equations," arXiv:gr-qc/0302071. 
[11] S. Frittelli and R. Gomez, "Boundary conditions for hyperbolic formulations of the Einstein equations," arXiv:gr- 
qc/0302032. 

[12] J.M. Bardeen and L.T. Buchman, Phys. Rev. D 65, 064037 (2002). 

[13] B. Gustafsson, H. Kreiss, and J. Oliger, "Time dependent problems and difference methods," John Wiley & Sons, New 
York (1995). 

[14] G. Calabrese, J. Pullin, O. Sarbach and M. Tiglio, Phys. Rev. D 66, 041501 (2002). 

[15] P.D. Lax, and R.S. Phillips, Commun. Pure Appl. Math. 13, 427 (1960). 

[16] P. Secchi, Diff. Int. Eq. 9, 671 (1996); Arch. Rat. Mech. Anal. 134, 595 (1996). 

[17] J. Rauch, Trans. Am. Math. Soc. 291, 167 (1985). 

[18] H. Kreiss, Commun. Pure Appl. Math. 23, 277 (1970). 

[19] A. Majda and S. Osher, Commun. Pure Appl. Math. 28, 607 (1975). 

[20] L.E. Kidder, MA. Scheel, and SA. Teukolsky, Phys. Rev. D 64, 064017 (2001). 

[21] A. Anderson and J.W. York, Jr., Phys. Rev. Lett. 82, 4384 (1999). 

[22] S. Frittelli and O.A. Reula, Phys. Rev. Lett. 76, 4667 (1996). 

[23] H.O. Kreiss, J. Lorenz, "Initial- Boundary Value Problems and the Navier-Stokes Equations," Academic Press, (1989). 

[24] When discussing initial-boundary value problems one sometimes demands a slightly stronger estimate that also bounds 
the norm of the solution at the boundary surface. For the purpose of the present article, it is sufficient to consider the 
weaker estimate (4) since we will show that when the determinant condition is violated, the inequality (4) cannot hold for 
all initial data 

[25] This condition is weaker than the uniform Kreiss condition [18] that requires that | det L_ | must be bounded away from 
zero. The reason why here we do not require the uniform Kreiss condition is that it might be too strong for the case of 
CPBC in General Relativity. As we will see in section IV the well posed CPBC that were derived in Ref. [7] do not satisfy 
the uniform Kreiss condition. 

[26] Notice that the characteristic variables defined here are related with the ones defined in Ref. [7] according to — 

[27] For n = 2 or n = 8/3 these fields are not complete. When r\ — 2 it follows that dj is traceless and as a consequence, 
8 AB Vfl = 0. But in this case one has the additional field 2C + C xx that propagates with zero speed. When n = 8/3, 
Vij — can be replaced by the fields Vab and 6V x a + 5C x a- 

[28] Actually this procedure is not unique since there is an ambiguity when first derivatives of the variables Kij and fkij are 
substituted for second derivatives of the three-metric. This ambiguity stems from the fact that one can always change 
the resulting expression by using the constraints Cikij = 0. For definiteness, we take the choice that leads to the same 
boundary conditions as in Ref. [7] . 



